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Abstract. 

We consider a few cases of homogeneous and isotropic turbulence differing by the mech- 
anisms of turbulence generation. The advective terms in the Navier-Stokes and Burgers 
equations are similar. It is proposed that the longitudinal structure functions S'„(r) in ho- 
mogeneous and isotropic three- dimensional turbulence are governed by a one- dimensional 
equation of motion, resembling the ID-Burgers equation, with the strongly non-local pressure 
contributions accounted for by galilean-invariance-breaking terms. The resulting equations, 
not involving parameters taken from experimental data, give both scaling exponents and am- 
plitudes of the structure functions in an excellent agreement with experimental data. The 
derived probability density function P{Au,r) ^ P{—Au,r) but P{Au,r) = P(— Am, — r), 
in accord with the symmetry properties of the Navier-Stokes equations. With decrease of 
the displacement r, the probability density, which cannot be represented in a scale- invariant 
form, shows smooth variation from the gaussian at the large scales to close-to-exponential 
function, thus demonstrating onset of small-scale intermittency. It is shown that account- 
ing for the sub-dominant contributions to the structure functions Sn{r) oc r^" is crucial for 
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derivation of the amplitudes of the moments of velocity difference. 
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Introduction 

Intermittency of turbulence, not contained in the Kolmogorov theory, is one of the most 
intriguing and mysterious phenomena of continuum mechanics. Experimentally detected in 
the early sixties, this feature of high Reynolds number turbulent flows still remains a major 
challenge to turbulence theory. Landau's 1942 remark [1] that the large- scale fluctuations of 
turbulence production in the energy-containing range can invalidate the Kolmogorov theory 
was one of the motivations for construction of various "cascade" models attempting to explain 
this phenomenon, manifested in the anomalous scaling of the structure functions 5'„(r) =< 
{u{x + r) — u{x))'^ >=< >oc AnV^" with the exponents ^„ 7^ n/3. The first model 
of this kind was proposed by Kolmogorov himself in 1962 [2]. Recently, some important 
analytic advances, leading to evaluation of both scaling exponents and the amplitudes An 
were made for the problems of the passive scalar, advected by a random velocity field and 
the random-force-driven Burgers turbulence [3] -[7]. First, it was proposed by Kraichnan [3] 
that scalar structure functions < {T{x) — T{x -|- r))^" > can be solutions of homogeneous 
differential equations , thus leading to the non-trivial values of the exponents ^„ which 
could not be found from dimensional considerations. Then, Gawcdzki and Kupiainen [4], 
Chertkov et. al. [5]-[6] and Shraiman and Siggia [7] showed that, indeed, it was the zero 
modes which were responsible for the anomalous scaling in some limiting cases of the passive 
scalar problem. Similar results were arrived at in many of the following studies [8]- [10]. 

Polyakov's theory of the large-scale-random-force-driven Burgers turbulence [11] was 
based on the assumption that weak small-scale velocity fluctuations {\u{x+r) —u{x)\ << u„ns 
and r « L, where L is the the energy input- scale), obey the galilean invariant dynamic 
equations, meaning that the integral scale and the random-force- induced single-point Urms 
cannot enter the resulting expression for the probability density having the scale-invariant 
form: 

P([/,r) « -F(-) (1) 
r r 

When U > Urms galilean invariance (GI), even of the small-scale dynamics can be violated 
and the PDF scales with U /urmsi leading to saturation of the scaling exponents ^„ = 1 forn > 
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He = 1- This general feature of Burgers turbulence was confirmed by numerical experiments 
in which turbulence was generated by different random forces generating various scaling 
exponents ^n<nc [12]-[14]. It was shown in Ref. [12]-[14] that the value of the critical moment 
number ric depended on the forcing function spectrum. Still, at n > all Sn oar indicating 
that decorrelation introduced by the noise was too weak to prevent the shock formation. 
Recently Chertkov ct. al. [15] obtained a similar structure of the theory considering a one- 
dimensional problem of a passive scalar advected by a random velocity field. In the limit 



where, to simplify notation, we have set the values of all numerical constants equal to unity. 
Similar result was also obtained in the work on the large-scale- driven Burgers turbulence in 

when space dimensionality — oo by Bouchuod et. al. [16]. 

Polyakov's idea [11] about the role of violation of galilean invariance in generation of anoma- 
lous scaling resonates with Landau's remark about important influence of the large-scale 
fluctuations on the small-scale dynamics [1]. In this paper we will attempt to combine the 
zero-mode and the Gl-breakdown ideologies to derive equations governing the probability 
density function of the longitudinal velocity differences in strong turbulence. We would like 
to reiterate that two measurements in the same turbulent flow performed in the laboratory 
and in the moving frame of reference (train or ship) must give the same answers. The "vi- 
olation of galilean invariance" is understood hereafter only in a limited Polyakov's sense as 
a possibility of Urms entering the probability density of velocity difference. It will be shown 
that the details of the large-scale turbulence production mechanism are important, leading 
to the non-universality of probability density function (PDF) of velocity difference. The 
results will be compared with experimental data. 
Formulation of the problem 

Turbulence in Nature results from hydrodynamic instabilities of various laminar low- Reynolds- 
number flows. The transition phenomena are is not universal, depending on geometry, ex- 
ternal fields etc. and, at the present time, cannot be accounted for by turbulence theory. 
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Hoping for some universality of the small-scale velocity fluctuations in the inertial range, it 
is customary to develop a theory of turbulence, driven at the large scales by some terms in 
the Navier-Stokes equations, one can treat theoretically. Usually, these large-scale forcing 
terms are assumed to be irrelevant in the inertial range. Below we discuss three models 
corresponding to different mechanisms of turbulence generation. First, let us consider the 
Navier-Stokes equations on an infinite domain: 

Vt V • Vv = f - Vp i/VV (2) 
V • V - 

The gaussian large-scale forcing f is defined by the two-point correlation function 

< /,(x, t)/,(x', t') >= A,«;(|x - x'|)5(t - t') (3) 

where the projection operator Pij ensures the incompressibility of the solution. The force is 
assumed acting in the interval of wave-numbers << ko ~ 1/L >> kd where Vt oo 

and 77 = l/kd are the volume of the system and dissipation scale, respectively. In other 
words, the forcing spectrum is assumed to decrease very rapidly outside the interval k ^ ko- 
In the limit k « ko the system is in thermodynamic equilibrium and is described by the 
gaussian statistics and energy spectrum E(k) oc k'^ [17]. Thus, the order of the hmits is: 
first we set the large value of ko and then z/ 0, so that k^/ko — > 00. 

In the case of the 5-correlated forcing function (3) the source-related contribution to the 
equation for the two-point correlation function can be written in a very simple way: 

W =< Vi{x, t)fi{x\ t) >oc k{\x - x'\) 

2. Often, in the real- life experimental situations, when turbulence is generated in the 
vicinity of the boundaries (wall flows), nozzles (jets), bodies (wakes) and later on transported 
into the bulk of the flow where the measurements take place, the model (l)-(3) is not valid. 
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In this case a better description is given by the initial value problem, taking into account 
the turbulence decay during this delay time. This will be also discussed in what follows. 

3. In the third class of the flows turbulence is produced by a large-scale shear. Then, 
introducing the so called Reynolds decomposition 

U = V + V 

where v (< v >= 0) is the fluctuation from the time- independent mean velocity < U >= 
V, the equation for velocity fluctuations v is given by (1) with f = and the turbulence 
production term in the right side [1]: 

(v • V)V 

The force-free Navier-Stokes equations are invariant under rotations, space-time trans- 
lations, parity and scaling transformations. They are also invariant under galilean trans- 
formations X — > X -I- Vt and v — > v -|- V where V is the constant velocity vector of the 
moving frame. Boundary conditions and forcing can violate some or all of the symmetries 
of the equations (1). It is, however, usually assumed that in the high Reynolds number 
flows with 1/ — > all symmetries of the Navier-Stokes (even Euler) equations are restored 
in the limit r —> and r >> rj where rj is the dissipation scale where the viscous effects 
become important. This means, among other consequences, that in this limit the root-mean 
square velocity fluctuations Urms — V<v'^ >, not invariant under the constant shift, cannot 
enter the relations describing moments of velocity differences. If all this is correct, then the 
effective equations for the inertial-range velocity correlation functions must have the sym- 
metries of the original Euler equations. For many years this assumption was the basis of 
all turbulence theories. Based on the recent understanding of Burgers turbulence [11]- [14], 
some of the constraints on the allowed turbulence theories will be relaxed in what follows. 

We are interested in the multi-point velocity correlation functions: 

C„(xi,X2, ....X„) =< t'a(xi)t'j2(x2)....t'm(Xn) > 
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and longitudinal structure functions: 



5„(r) =< {{u{x + r) - u{x)y >=< (Aii)" > 

where u{x) is the a;-component of the three-dimensional velocity field and r is the displace- 
ment in the direction of the 

In 1941 Kolmogorov, considering decaying turbulence, derived an equation for 5*3 (r) valid 
when r ^ 0: 

ri dr r4 ^r^ ^ ^ 

leading to the famous Kolmogorov |-law: in the limit i/ — > and atL>>r>>77 where 77 
is the dissipation scale of turbulence defined as: 

vSrM = 0{l) 

The third-order structure function in the inertial range L>>r>>77is given by a 4/5 law:: 

o 

In fact, the expression (4) is an approximation, neglecting small contributions from time- 
derivatives and source functions, valid in the limit r — > 0. In general, it must be modified to 
include the forcing function: 



^3(r) = --^r + 0(r(AK)(A/)) (5) 
o 

We can see that for a white-in-time large- scale forcing function K{r) = «:(0) — 7r^, the 
sub- leading contribution to Kolmogorov 4/5-law in the inertial range r — > is 0{r^). In a 
more general situation evaluation of the correction is not so simple. 

In case of decaying turbulence the sub-leading contribution to (5) is : 

while when turbulence is produced by the large-scale shear, it is somewhat different [18]: 
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When r is small, these terms can be neglected. However, as will be shown below, they 
must be accounted for since the procedure of evaluation of the probability density P{Au,r) 
involves matching of the inertial range and large-scale solutions. 

Derivation of (4)- (6) is based on the fact that due to the incompressibihty condition, all 
transverse correlation functions can be expressed in terms of the longitudinal ones leading 
to the closed equations. One can say that in a very limited sense the procedure projects 
the original three-dimensional problem onto a one- dimensional one. Regretfully, due to the 
coupling between different components of the velocity field, caused by the pressure terms, we 
cannot rigorously derive similar expressions for the high-order moments Sn{r). The second 
difficulty is in the presence of the dissipation anomaly (see below) . Still, we can attempt to 
use some general features of the equations of motion and derive the scaling properties and 
general form of the probability density function P{Au,r). 
Equations for the probability density 
One can introduce a generating function 

where the vectors Xi define positions of the points denoted by the numbers 1 < i < N. Using 
the incompressibihty condition, the equation for Z can be formally written: 

— + -I +1 + D (7) 

where If,Ip and D are related to the forcing, pressure and dissipation contributions to the 
Navier-Stokes equation (see below). Although the advection contributions are accurately 
accounted for in this equation, it is not closed due to the pressure and dissipation terms. 
The latter can be treated using Polyakov's operator product expansion ideas [11], while the 
former presents an additional difficulty to be dealt with. In what follows we will be mainly 
interested in the moments of the two-point velocity differences which in homogeneous and 
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isotropic turbulence can depend only on the absolute values of two vectors (velocity difference 
v(x') — v(x) and displacement r = x' — x) and the angle between them with 6* = 7r/2 and 
^ = corresponding to transverse and longitudinal structure functions, respectively. In 
spherical coordinates the explicitly written advective terms in the equation (7) involve 

«(||)^ o^^y^ 

and various trigonometric functions and angular differentiations. The theory of the longitu- 
dinal structure functions, presented below, is based on the assumption, correct for the third 
order moment 5*3 (r) (see (4)), that the angular dependence can be accounted for in a simple 
way and, as a consequence, there exist an equation for ^ = 0. This assumption is supported 
by the following observations. It is easy to show that in the inertial range the second-order 
structure function 

with DiLir) —< {u{x) — u{x + r)y >. More involved relation can be written for the fourth- 
order moment [16]: 

S4r,e) = DLLLL{r)cos\e) - 3DLLNN{r)sin\2e) + DNNNN{r)sin\e) 

where Di^i^jsfN =< {v{x)—v{x+r))'^{u{x)—u{x+r))'^ > and v and w are the components of the 
velocity field perpendicular and parallel to the x-axis, respectively. As one can easily deduce 
from the angular dependence, the functions DiLiLij) and D^NNNiT) denote longitudinal 
and transverse structure functions, respectively. In the limit ^ — > 

S,ir,e)^Di^LLL{r)cos\9) + 0{e^) 

rapidly approaching S4{r,9 = 0) = Dllll^t) = 5*4 (r). Based on the above expressions, 
we conclude that, as in the theory of multidimensional Burgers turbulence [16], where the 
probability density of velocity difference can be represented as: P{U,r,9) ^ P{U,rcos{9)), 
here in the limit 6' — > the mixing of the longitudinal and transverse correlation functions is 
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very weak {0{9^)). As a consequence, we assume that the closed equation for the probabihty 
density of longitudinal velocity differences exist. Generalization of the theory to the case of 
an arbitrary (not small) angle 9 is the subject of an ongoing study. 

We select N points Xi with < i < N on the x-axis and introduce the longitudinal generating 
function for the A'^-point correlation function: 

Zn =< e^'"(^*) > (9) 

where = ikj. The equation of motion for can be formally derived (we neglect the 
subscript N in what follows): 



Substituting the N-S equations into this relation gives: 

Z.^Y: < e'^<^^^[-u{x,)^^ + Uxj)] >+I^ + I^ + D (10) 

■ axj 

where fx{xj) is the x-component of the forcing, xj is the coordinate of the j*'* point and ^ 
is the partial derivative in the x direction only. The summation in the expression (10) is 
over the positions Xj and over the Greek subscripts a — 1; 2 denoting the components of 
the velocity field in the directions perpendicular to the x-axis. The life-threatening terms in 
(10) are: 

Ir = E < e^-(-)[-^.(x,)^] > (11) 
= - E < e"'"(^')^^ > (12) 

and 

D = A,- < e^^"(-^)^^^ > (13) 
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The theoretical and numerical work [16] on the mult i- dimensional Burgers equation led 
to the probability density and moments of velocity difference basically independent on the 
space dimensionality: the moments Sn<i{r) oc r" while Sn>i oc r. This is an indication that 
the shock production, dominated by the longitudinal components of the nonlinearity UidiUi 
(no summation over the subscript i), prevails over the processes coming from the mixed 
terms of the kind UjdjUi which can be neglected. In other words the multi-dimensional 
Burgers equation is well approximated by the system of weakly interacting ID-equations 
acting along various coordinate axis. It is clear that geometry of the objects generated by 
this system is very complex. The recent paper by Gurarie and Migdal [16], dealing with 
the two- and three-dimensional Burgers turbulence, introduced an angle 6 between veloc- 
ity difference and displacement vectors v(x -|- r) — v(x) and r, respectively, and using the 
instanton formulation, derived an expression for the generating function (see below) Z2 in 
the form Z2 oc exp{{Xr)'^ f{cos9)) with 7 = 3/2 independently on space dimensionality and 
A = |A[ and r — |r|. The calculated function f{cos{6)) ensured correct angular dependence 
of multi-dimensional structure functions. When ^ = 0, the derived expression basically re- 
covered the one-dimensional Polyakov's result. This result tells us that it is the projection 
of the velocity field on the direction of the displacement vector r that produces dynamically 
significant contribution to the multi-dimensional structure functions and that the longitu- 
dinal structure functions in three-dimensional Burgers turbulence are close to those in the 
ID-turbulence. Indeed, numerical simulations [16] of the three-dimensional Burgers tur- 
bulence revealed very complex velocity field with the structure functions Sn{r) very close 
to the ones, previously obtained in one-dimensional simulations. The possible smallness of 
interaction between different components of the advective terms is only part of the story. 
The pressure contribution Ip leads to effective energy redistribution between components of 
the velocity field and plays an important part in the Navier-Stokes dynamics. The pressure 
effects are nonlocal, instantaneously transporting information between different, even very 
remote, parts of the fiow. That is why the effects coming from the boundary conditions 
and large scale forces cannot be neglected even in description of the small-scale phenomena. 
The possible spontaneous breakdown of galilean invariance is the central assumption of this 
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work. 

The equation (10) can be rewritten as: 



A2 7 1^7 

Z, = E A.A,.(|.. -.,\)Z-Y.\^^'-^^]^h^I, + D (14) 

where the large-scale gaussian random force is defined in the limit r ^ by the correlation 
function: K{r) — k{Q) — ^r^ . Approaching the integral scale L the force correlation function 
K{r) rapidly goes to zero. We will see that the equation for the probabihty density of the 
velocity difference P(f/, r) where U = u{x + r) — u{x) contains the combination /t(0) — K{r) 
which is large at the large scales r ^ L. That is why the large-scale dynamics, dominated 
by the forcing term, show close to gaussian behavior of at least the first few moments 

Evaluation of It, Ip and D in (14) is a difficult problem and we have to make some as- 
sumptions. It is seen from the definition of the generating function Zjq that when two points 
merge, i.e. Xi — > Xj the iV-point generating function becomes the N — 1-point generating 
function with A = Aj -|- A^. This means that if, for example, the equation for Z2 contains the 
two-point sum: 



(/?(Ai, w(xi),xi, -Q^)Z2 + </?(A2,xi(x2),X2, ^ a{X, y)(p{X, u{xi) , Xi, -Q^)Z\ 

Here ip depends on the structure of the equation of motion and y — \xi — Xj\ 0. We 
assume that in this limit the unknown function a(A, 0) is finite. Not all functions </? satisfy 
this equation. For example: ip — \ and </? = A^, do. The functions ip can also include space 
derivatives. Combined with the general symmetry properties of (14) we can narrow the class 
of possible solutions and derive equation for the probability density. It follows from the 
Navier-Stokes equations that the theory must be invariant under transformation: A —A 
and Xj — > —Xj. In what follows we will adopt Polyakov's result that the main effect of the 
longitudinal part of the dissipation term is a renormalization of the coefficient in front of 
the 0(1/Aj) terms in the right side of (14). Based on the above considerations, we can write 
two models for the A^-point generation function, corresponding to the different mechanisms 
of turbulence production introduced above (see below) : 
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ff'7 h 87 8 
= X X.^,m.. - .,l)Z - Elgp^ - + C i: A,^Z + /. + + C (15) 

for the cases 2 and 3 and 

Z. = E A.A,.(|.. - - Elap^ - J-g^ ^Ir^I^^D' (16) 

J J J J J 

to describe turbulence generated by the white-in-time forcing (3). Here D' involves only 
"transverse" components of the D-term defined below. The 0(A^)- term in (15) comes 
from the 0(v) turbulence production, while K{r) — K{0) — /3r^, leading to the negligibly 
small O(r^) sub-leading contributions to the velocity difference structure functions, ensures 
the close-to-gaussian single-point probability density. The equations (15)-(16) violate neither 
"fusion rules" introduced above nor general symmetries of the Navier-Stokes equations. By 
dimensionality, the coefficient C/k{0) — 0{l/ul^g). The 0{K{r)) term in the equation (16), 
stems from the forcing function (3). Below we will discuss the two cases in detail. 

If one is interested in a single-point probability density, the equation (15) is to be solved 
for all Xj = X and A = Y^j^j- All space derivatives disappear due to homogeneity of 
turbulence and we have an expected result [17]: 




This fixes the value of the coefficient C . Thus, the equation (15) yields gaussian distribution 
of the single-point velocity field. This result will be used below as a matching condition for 
the probability density P{U,r) in the large-scale hmit r ^ L. The experimentally observed 
single-point probability density P{u) is very close to but not equal to gaussian deviating 
from it at the large values of velocity fluctuations u » Urms- The theory, developed here, 
is applicable to any expression for P{u), not only to the gaussian. Still, the gaussian, which 
will be used below to compare the theoretical predictions with the data, is a very good 
approximation. 

We need an equation for the the generating function Z2 with Ai -|- A2 = giving: 
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Z2{X,r) =< e^W^+'-)-"(^)) =< e^^ > (17) 
In a statistically steady state: 

- = ''(^^■(°' - ^^'(l^- - + It + I, + D-- (18) 

Where 



/. =< A[^ - Ml^ii) ]eA(«(2)-(i)) A < ve^^ > 
oxi 



and 



i?' = A^ < [9^,„i.(x2) - 9^^„^(xi)]e^("(^)-«(^)) >^ A < d'e^^ > (19) 

The point merging in It and Ip must be regular while the same procedure in the dissipa- 
tion term D' involves divergences which are cancelled by viscosity i/ — > 0. As was pointed out 
in [11] the longitudinal, 0(9^ii(x)), components of the D-term result in the renormalization 
of the coefficient in front of the 0(1/ A) - contribution to the right side of (14). We are still 
left with the remaining 0{d'^u{x.)) - piece of the dissipation anomaly, pressure terms Ip and 
the It - contributions, mixing all components of the velocity field. Having in mind general 
fusion rules, considered above, and the fact that the equation is invariant under transforma- 
tion A — > —A and r — > — r we, not being concerned with preservation of galilean invariance, 
write the equation for Z2 corresponding to (15): 

dXdr X dr r dX L dX ^ ^ 

The equation (20) includes the derived above expression for the coefficient C and the un- 
known parameters B° and A to be determined from the theory. The characteristic time in 
(20) T ^ L/urrns = 0(1) is independent on the displacement r. A natural generalization of 
this model is (20) with the last term in the right side: 
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T(r) dX 

with T(r) oc r^'^ with the exponent depending on the physics of the problem. In Kolo- 

mogorov turbulence C,t ~ 2/3. 

The model corresponding to (16) is: 

^_B^dZ2^AdZ2^ ^2^2^ .21) 

dXdr X dr r dX 
where the 0(r^)- contribution is to be kept. The equations (20)-(21) are based on the assump- 
tion that the dynamic role of the pressure and dissipation terms is in the renormalization of 
the coefficients in front of the already- present advective contributions (8) to the equation 
(7). Similar assumption was fruitful in the theory of Burgers turbulence [11]. Except for the 
last terms in the right side, these equations are the same as the one in Polyakov's theory of 
Burgers turbulence with the B° - term simply renormalizing part of the advection effects, 
present in the original equations. The meaning of the A-term will be discussed in detail 
below. It will be shown that it is not only responsible for prevention of the shock formation 
but it also makes the weak {\U\ « Urms ) structures of the Navier-Stokes dynamics much 
stronger than their counterparts of Burgers turbulence. It is clear that due to homogeneity of 
turbulence, all space derivatives and related to them ^-contributions go to zero when r > L 
. In the one-dimensional case the dissipation anomaly D, discussed in detail by Polyakov, 
leads only to a relatively small modification of a corresponding coefficient . It will be shown 
that in case of three-dimensional turbulence — 20, part of which is to be attributed to 

a large pressure effects preventing the shock formation. In the resulting dynamic equations 
(20)-(21) the contributions from D and Ip are mixed and the origins of each term, hidden in 
numerical values of the coefficients 6, B° and A, are not easy to establish. The equation (20), 
explicitly involving the single-point Urms is suited for description of the generating function 
in both r — > and r ^ L limits. In the inertial range, where the displacement r is small, the 
0{urms) and the forcing contributions can be neglected. They are important, providing the 
large-scale gaussian matching constraint, needed for determination of the amplitudes of the 
structure functions Sn- Thus, the equations (20)- (21), describing the correlation functions 
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in the inertial range, differ by the last in the right side terms reflecting the details of the 
large-scale turbulence generation processes. It will be shown below that this difference is 
responsible for non-universality of the probability density function of velocity difference. 

In the limit r — > the equation for the probabihty density is derived from (20)- (21) 
readily: 

d dP dP Ad Urms d'^ , . 

d ^BP ^BP A d ^ ^d^Z 
Properties of the solution 

Multiplying (22) by C/" and assuming existence of all moments leads to: 

dSn An Sn , Urms n^i'lT' ~ ^) ct I \ inA\ 

^ = + -r ^24) 

where B — —B° > 0. The equation (24) is to be solved under the constraint (4) which is 
the result of the energy conservation inherent to the Navier-Stokes equations. This is the 
consequence of the renormalization ideology leading to (24), which is a model, not rigorously 
derived from (l)-(3), but based on some general symmetries of the Navier-Stokes equations. 
That is why the 4/5 - law comes out of (24) only for a particular set of parameters. In a 
"final" theory the rigorous equation for 5'„(r) must automatically produce the correct result 
for S3{r). Neglecting the last term in the right side of (24) (see below) the solution for S3{r) 
in the limit r — > is derived readily: 

3A 

This means that the coefficient A = and = — |^^- Seeking the solution in the form 
Sn oc r^" we obtain: 



{3 + B)n 
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The normal Kolmogorov scaling = n/3 , corresponding to the no - intermittency case is 
achieved in the limit B ^ oo. The maximum, Burgers- like, intermittency with all exponents 

— i due to the tanhx-shocks is recovered when B — 0. 
Deriving (25) the contribution of the order 

L n + B 

was neglected in comparison with the 0{ ^^^^y ) terms. Substituting the expression for 
Sn — Anr^"^ with ^„ from (22) gives a general solution: 

5„=^„.^-+!^^„_. "'"-j) (26) 

L n + B t,n-l - ?n + 1 

The coefficients will be derived below. The expression (26) shows that due to the 

presence of the 0{r^"~'^~^^) subdominant contributions, the experimental determination of the 

scaling exponents is a difficult task and that proper accounting for it can lead to substantial 

broadening of the inertial range and more accurate determination of the numerical values of 

the scaling exponents In addition, it establishes the relation between the amplitudes of 

the odd and even-order moments. 

By definition of the integral scale, adopted in this work, the third-order moment 6*3 (r = 

L) = 0. Since Urms ~ i^L)^, the expressions (25)- (26) give: 

4 B + 2 

^3 = -- = -9A2- 



5 "(5 + 3)2 

Taking in accord with the closures derived from various renormalized perturbation expansions 
A2 ~ 2 leads to S fa 20 and: 

e. = (27) 

^ 3 n + 20 ^ ^ 

The calculated values of the exponents ^1/5 = 0.0759; ^1/4 = 0.0946; ^1/2 = 0.187; ^1 = 
0.365; ^2 = 0.696; ^4 = 1-278; ^5 = 1-533; = 1.769; ^7 = 1-988; = 2.190; ^9 = 
2.379; ^10 — 2.555 are indistinguishable from the best available experimental data. One has 
to keep in mind that the value of the parameter B can be non-universal slightly varying from 
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flow to flow. This can lead to some non-universality of the exponents. The comparison of 
the magnitudes of the exponents, given by (27), with the outcome of numerical simulations 
by Chen [19] is presented on Fig.l. The expression (27) predicts saturation of the values 
of the exponents ^„ at ^oo ~ 20/3 in agreement with some general ideas based on the path 
integral representation of the solution of the passive scalar problem Chertkov [20] . 
The expression for Sn{r) corresponding to the model (23) is: 

, , . n(n-l)(n-2) r^+^ns 
Sn = A^r^^ + 7A„_3^ ^ — (28 

' n + B 3 + -Cn ^ ^ 

It follows from (28) that 



in accord with an exact result. 

To calculate the value of parameter B one has to solve the equation for the probability 
density P{U, r) > going to gaussian in the hmit r/L — > 1. The equation is: 

^U?f - = k'-^)^U P (29) 

dU dr dr 2, ' dU ^ ' 

It may be somewhat easier to deal with the equation for the generating function Z2: 

a^Za _ B_dZ2 _ I ?> + B dZ2 
dXdr X dr r 3 dX 

The structure of the solution is clear from the scaling of the moments derived above: 

°° (3+B)n 



with as yet unknown amplitudes An > which will be evaluated below. The most important 
outcome of this expression is the fact that the odd-order moments S2n+i < which means 
that the PDF P{U, r) ^ P{-U, r). It is clear that P{U, r) = P(-C/, -r) in accord with the 
symmetry of the Navier-Stokes equations. 

To evaluate the probability density function we need to match the inertial range PDF 
with either the energy containing or dissipation range probability density functions. Based 



on the result for the single-point PDF [17], one has to seek the solution to (25) that becomes 
very close to gaussian at the scales larger than some integral scale L . This condition can 
serve as a definition of the integral scale. In reality, the integral scale is never much smaller 
than the size of the system. The experimental data show that at the large scales the PDF is 
close to gaussian with some small deviations seen at the far tails where U » Urms- Based 
on the data we can safely assume that at the large scales the firstt few (10-20) moments are 
very close to their gaussian values. The non-zero value of the odd-order moments 5'2n+i(r) 
with n > 1 implies the asymmetry of P{U,r). However, this asymmetry is very small with 
S2n+i{r) / S2n+i{r) « 1 where S2n+i(r) =< \u{x) — u{x -|-r)p"+-'^ > is often measured by 
the experimentalists. It has been shown [21] that up to n ~ 5 this ratio is in the range of 
0.1 and that the experimentally observed PDF can be made symmetric by a signal-filtering 
procedure, leading to near- vanishing of the odd-order moments. The procedure left the even- 
order moments unchanged, indicating that the PDF asymmetry only weakly infiuences the 
even-order moments. The data will be presented below. 

Now we set L = 1 and, neglecting all odd-order moments, calculate the PDF directly 
from the generating function: 

2 /■°° 

Pe(U,r) = - cosikU)Mk,r) (30) 

TT JO 

with 

oo jL2n 

0(A;, r) = 1 + $:(-l)M^(2n - 1)!! r«»— (31) 

n=l 

We see that this expression gives S2n = ^2(2^ — 1)" f^^" leading to the desired gaussian 
values at r = L = 1. The same result P(U,r) oc exp(—-^) is recovered in the limit B ^ 00 

^ ' ^ 2r3 

. In the opposite limit S — > the probability density tends to 



9n^ 



PeiU.r) = {l-r)6{U) + r^-e 

giving all ^„ = 1. This corresponds to Burgers turbulence in the Gl-broken range [11]. This 
result is very close to the outcome of the theory of the Burgers turbulence in the limit of 
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space dimensionality D ^ oo [16], predicting P{U,r) = (1 — r)5{U — r) + ripi^^^). This 
fact is an indication that the gahlean-invariance-breaking terms in the equations of motion, 
obtained in this work, are quite close to the truth. It is clear that to reproduce the shifted 
5- function of Ref.[16] we have to abandon the simphfication of treating the PDF as an even 
function of U. 

To uncover the inner structure of the (5-function the data are to be presented in coordi- 
nates [//r" with a ~ 0.365 (see Figures below). We can compare the prediction given by 
(31) with the experimental data on 

C2n{r) = = (2n - 1)!! r«--"«^ (32) 

where all ^„ are given by (22). Sreenivasan et. al [21] measured C2n(r) in the low Reynolds 
number experiment {Rx 200) conducted in the laboratory boundary layer. The results of 
Ref. [21] for r 0.1 - 0.2 are: 

C4f«3.6; C6~24.8; Cg ~ 261; Cio 3770 
which are to be compared with our predictions r = 0.2: 

C4 = 3.16; C6 = 25.05; Cg = 272; Cio = 3256 : 

The inter mittency grows strongly with decrease of the displacement r. For example: for 
r = 0.1 the relations derived in the present paper give: 

C4 = 3.42; C6 = 31.26; Cg = 412; Cio = 6184; 

These numbers agree extremely well with the results of numerical simulations by Chen 
[22] who was able to produce the data set at Rx = 200 consisting of a few billion points. 
The somewhat lower values of Cio obtained in the physical experiments can be attributed to 
the insufficient statistics of the real-life data set. Recent high Reynolds number experiments 
{Rx ~ 15000) [23] produced similar results for the moments Sn with n < 6. The comparison 
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between theory and experiment is somewhat difficult due to an uncertainty in the theoreti- 
cally needed ratio r/L where L is a poorly defined integral scale of turbulence. The inertial 
range prediction prediction of this paper C4 ~ 3r"°'^^^ agrees well with the data obtained in 
the high Reynolds number experiments {Rx ^ 1500 — 2500) in the planetary boundary layer 
[21]. 

The evaluated probability density functions are compared with the outcome of the mea- 
surements of NouUez et. al. [24] on the Figs. 2-4. The high quality experimental data on the 
transverse structure functions were obtained using real-space measurements in the air jet. 
Thus, for a time being, comparing theory with experiment here we assume that the proba- 
bility density of the the longitudinal velocity differences has the same shape as the one of 
transverse velocity differences. The quantitative agreement for 1/4 < r/L < 1 is very good. 
We were not able to plot P{U,r) for very small values of the displacement r but the data 
on Fig.4 show the tendency of the PDF to the (5-function in the limit r — > in accord with 
the analytic asymptotics. The Fig. 2 presents the same curves, plotted in the coordinates 
U/r^/^. One can see increasing deviations from the gaussian at r/L = 1 with decrease of 
the displacement r. It is clear from the Figure that the probability density cannot be rep- 
resented in the scale-invariant form (1). This is the manifestation of the intermittency. The 
calculation was performed using Mathematical'^ which failed to produce the generation 
function (f){k, r) with k > 5. The information on (f){k, r) with < A; < 5 was sufficient to 
calculate P{U, r) with the accuracy fa 1 — 2% in the range of variation < [/ < 3 — 4 and 
1/4 < r < 1. The fact that at r = 1 the PDF must be gaussian serves as a good test of the 
quality of the numerical procedure. 

Using the derived expression for Z2, one can easily evaluate correct asymmetric probability 
density P{U,r) giving the values of the odd-order moments in agreement with experimental 
data. All one has to do is to introduce a term involving the szn-contribution to the Fourier- 
transform in (26)- (27), generating non-zero odd-order moments. Demanding that S2n+i{L) ~ 
we have from (26): 

^ _ Urms 2n{2n + l)A2n 

~ {£L) 3 (2n + 1 + 5) (e2n - 6n+l + 1) 
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valid for n > 1. Setting Urms ~ (^-^)^ ^^'^ taking /I2 ~ 2 we obtain: 



^3 ~ -0.8; A5 ^ 23.0; A7 ^ 650 : 

In the case of the large-scale- driven turbulence the expression for A2n+i is derived in a 
similar manner: 



_ 2n(2n+ l)(2n- 1) 



'(2n+l + B)(3 + 6n-2-5n+l) 

This relation contains two unknowns 7 and B. Assuming universality of the exponents, 
implying B — 20, we find from the relation for A3 = —4/5 that 7 ~ 6. In the correct 
dimensional units 7 ~ 6£^. This will be important below for the quantitative comparison 
of theoretical predictions with experimental data. This result has some interesting conse- 
quences. Keeping the amplitudes equal to the ones derived above: 

A5 14; A7 fa 373; Ag ^ 28500; 

Wc see that even if the exponents arc universal, the amplitudes of the moments are not, 
meaning that the shape of the probability density can vary from flow to flow. It is interesting 
that to experimentally observe this effect one has to measure high-order moments since the 
first few structure functions in different flows, considered here, seem to be close. Hereafter 
we will be mainly interested in the model (20) which is more closely related to physical 
experiments. 

The expression 

is to be substituted into the integral: 



o 

TT 



/oo 
sin{kU)(f)i{k,r) dk (34) 
-00 



to give a total asymmetric PDF P{U, r) = Pe{U, r) + Po{U, r). A very accurate parametriza- 
tion of the central part of the PDF, corresponding to not-too-large values of U, illustrating 
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appearence of the asymmetric PDF is: 



Mk,r) ^ -^rk^(l){k,r) (35) 
io 

giving: 

2r 8^ P (U r) 

P{U, r) = Pe(C/, r) + P,(C/, r) = Pe(C/, r) + - ^^3' ^ (36) 

This expression is approximate and serves only to illustrate the mechanism of appearence of 
the experimentally observed asymmetric PDF P{U,r). The central part of the probability 
density can can be very accurately parametrized by the formula (36) with: 

Pe(t/,r) = ^e--*^°'^? (37) 

where N is a, normalization constant, x ~ ^ with, as above, a ~ 0.365. The parameter 
a 2 — 4. The expression (31)-(32) gives a very good approximation for the moments 5'„(r) 
with n < 10. For example: A2 = 1.9-2.6; A3 = -0.8; ^4 = 20-28; A5 = -15.4-21; Aq ^ 
657; At —785; etc., very close to the first few amplitudes An, calculated above. Figs. 
5-6 show the probability density function P{U,r) given by (36)-(37). 

The theory can be approximately generalized to the case of correct anomalous exponents: 

-Rn = — ~ -^r'^ 2 (38) 

A2 

The expression (38) gives 5*3 in accord with the Kolmogorov relation. For r f=:i 0.1 and A2 ^ 2 
we derive: P3 = -0.31; R5 = -4.2;. 

To assess internal consistency of the theory and understand the role and meaning of 
the ^-contribution to (20), we write a formally exact equation for the two-point generation 
function: 



d^Z2 2 dZ2 



X\k{0) - K{\xi - Xj\))Z2 + IT + Ip + D (39) 



dXdr X dr 

The equation for the probability density P{U,r) can also be formally written: 
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where 



d dP dP 1 9^ 



(40) 



T(t/,r) =<ip + iT + d\U > 



(41) 



is the conditional expectation value of T = ip-\-iT-\-d for a fixed value of velocity difference U. 
Here d — d' -\-v{ux'x'{x') —Uxx{,x)). The formulation of probability density functions in terms 
of conditional dissipation and production was introduced in [25]- [26] and became a subject 
of intense investigations. Three-dimensional turbulence problem is somewhat different since 
there is no way one can separate various contributions to T{U,r). Indeed, we are dealing 
here with the projection of a three-dimensional dynamics onto a line. This leads to violation 
of all conservation laws and it is clear that the pressure-induced processes of the correlation- 
redistribution between different components of velocity, probably conservative in the original 
equation of motion, can lead to the dissipation-like effects of T{U,r). We can see from the 
above definitions that due to homogeneity of the turbulence 



It will be shown below that (42) can serve as an equation for determination of the only 
unknown parameter of the theory B. Comparing (40)-(41) with (22) gives: 



-l^TiU,r)P{U,r)^-^^UP{U,r) + {B'^-l) f_J-^^dy-A^P{U,r) (43) 



RecaUing that B = -B° > and ^ = (3 + 5)/3 > 0, the relations (42)-(43) give an 
equation for the only unknown coefficient i? > 0. It can be solved numerically to establish 
consistence with the above calculation leading to S ~ 20. An interesting insight into (42)- 
(43) is obtained if we use the fact that the deviations from scaling of a central part of the 
PDF are small. This means that: 



iT-\-ip-\- d — 



J —oo 



/ T{U,r)P{U,r)dU = 



(42) 



P{U, r) 
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Substituting this into (42)- (43) gives: 



Now, we have to define the "representative exponent" a ~ ^n/n with < n < 3, which can 
be done only approximately. Choosing a — and recalling that ii{B + 1) = ^4, the last two 
terms cancel and 

T{U,r) = 2^U 
ij 

satisfying the constraint (42). 

A much more interesting relation is obtained by substituting (36) into (43) with the 
scale-invariant expression for Pf.{U,r). Taking into account that a{B — A, we have: 

nU.r) ^ 2(g + !)(!- 3a) dPe{U,r) 

2 15P(C/,r) dU ^ ' 

The plot of T{U, r) calculated from (39) using an approximate relation (36)- (37) is presented 
on Fig. 7 for r « 0.01. Outside the interval —2 < f//r" < 3 — 4 the curve saturates 
which might be an artifact of the exponential asymptotics the approximate formula (36)- 
(37), valid only when U not too large. It is interesting that the curve is asymmetric, reflecting 
the asymmetry of the PDF. In the limit C/ — > the expression (44) gives: 



15 r2" 

This relation takes into account that 3q; > 1. Choosing a || 0.383 gives: 

m.)» 0.84-^. (45, 

where ^2 ~ 2 is the amplitude of S2{r). 

Coming back to the estimate of parameter B. We saw that setting a — ^i, though 
shedding some light on the structure of the expression (43), does not allow it's estimate. 
The problem is that Si = and one cannot use it for the non-dimensionalization of the 
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argument of the scale-invariant PDF. A different parametrization a = .^2/2, typically used 
for analysis of experimental data, leads to cancellation of the last two terms if 



«i(B + l) = £±l 

2 ^ ^ 3 

giving B 40. With the exponent a ~ 0.383, best characterizing the top of the PDF 
Pe{U, r) we have B ^12. This simple calculation demonstrates consistency of the model for 
T{U) with the magnitude of the parameter B 20, derived above. 
The expression (43) can be modified using an identity [27]: 

^^U^r) ^ --i-r<^\U> P{U,r) ^ i-h{U,r)P{U,r) 



dr dU dr dU 

dU 
dr 



where Ii{U, r) is the conditional expectation value of for a fixed value of velocity difference 



U. The equation (43) reads: 



^T{U,r)P{U,r) = -^-^UP{U,r) + {B + l)h{U,r)P{U,r) - A^P{U,r) 



It is easy to see that 

IdU u{r + 6) - u{6) - u{r) + u{0) 



^u,{r)-u,{0) (46) 



2dr S 
This expression shows that experimental and numerical investigations of conditional expec- 
tation value of velocity-derivative difference for a fixed value of U is important since the 
combination 

[(B + l)<f\U>-^'^]P{U,r) (47) 

determines the structure of the probability density. Since P{U,r) — P{—U,—r), we have 
from (35): 

T{U, r) = -T{-U, -r) 
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Large-scale corrections to scaling 

To complete comparison of the present work predictions with experimental data let us discuss 
some consequences of the expression (26). Experimental determination of the inertial range 
is usually done by establishing the interval where the third-order moment S^lr) oc r, in 
accord with the Kolmogorov relation. As one can see from (26) it is not so easy because of 
the 0{rS2) subdominant contribution giving 

S3 ^ -0.8r + 0(r^-^) 

where we set ^2 ~ 0.7 and Urms = L = S = 1. To make quantitative comparison of this 
relation with the data we have to restore physical dimensional units. Let us define 

Snir) = AnSn = ^„(£:r) 3 (^)^"- 1 
The relation (26) can be rewritten: 

The high-Reynolds number experiment [23] was conducted in the atmospheric boundary 
layer on the tower at 35m above the ground. The measured Urms ~ 1 Am/ sec and the mean 

2 

dissipation rate S ~ 0.03^3. The measured root-mean-square stream- wise velocity compo- 
nent Urms in the wall-bounded flows is somewhat larger than that of the velocity components 
perpendicular to the direction of the flow. Since in the relation (26) we are interested in 
Urms: Corresponding to the top of the inertial range, a good estimate for Urms/{^L)^ fti 1. 
Taking L Ri 35m and A2 ~ 2.0 — 2.5 gives: 

c 

--^ fa 0.8 - (0.06 - 0.08)r°-^ 
tr 

and 

_:?^~20 - 2r'^-^5 
The experimentally observed [23] relation 
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,^5 = ^(£r)i(y)«^-i ^ 0.09r^-5^ 

Ij 

is consistent with numerical values for £ and L used above. The third and fifth-order mo- 
ments, calculated from the above relation, are presented on Fig.8. The parameters used were: 
Ki 20 and A4 ^ 22 [23]. The value of the integral scale L a; 35m, used above, can be a bit 
overestimated. Choosing L 20 — 30 m does not substantially modify the above conclu- 
sions. This result is extremely important since it shows that without explicit accounting for 
the subdominant 0{r^''^)- component of S3, one cannot observe the Kolmogorov relation but 
in very high Reynolds number flows. It also tells us that, in fact, inertial range can be made 
much broader and the scaling exponents can be established very accurately with the proper 
data processing. One can see from the Fig.8 that the fifth-order moments starts deviating 
from its asymptotic value earlier than ^'3. The expressions corresponding to the model (23), 
can be written easily: 

^^n+i _ . 7A2„-2(2n + l)2n(2n - 1) r .3+6n-2-6n+i 

S2n+i n + B L 

As we see, in the limit r — > the contribution from the sub-leading term in this case is 
much smaller. This means that for a given Re investigation of the scaling exponents ^„ in 
the numerically- created large-scale-driven turbulence is much easier. 
Discussion and Conclusions 



The theory developed in this work is based on equations including a simple model for 
the pressure and dissipation terms. This model, though satisfying all basic symmetry con- 
straints, has not been rigorously derived from the Navier-Stokes equations. The merits of 
such work can be judged by comparison of theoretical predictions with experimental data. 
The calculated exponents and the amplitudes of the structure functions 5'„(r) agree very 
well with available experimental data. The theory also predicts the large-scale corrections 
to scaling, thus allowing calculation of S'„ up to the large-scale cut- off L at which S'2„+i(r) 
become very small. This prediction is non-trivial and can serve as a rigorous test of the 
model. The scale where it happens is an integral scale of turbulence, corresponding to the 
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top of the inertial range. This definition seems very plausable since it corresponds to the 
length-scale of the non-zero energy flux set up. 

The most straitforward experimental test of this theory can be performed in a following 
way. Assuming that the odd-order moments have the form: 

S2n+1 = A2n+ir^'"+' + B2n+ir'''-+' 

where the exponents reflecting the large-scale dynamics, have been avaluated above for 
the two cases of turbulence production. The Log-Log plotting of the functions 

^2n+l - B2n+ir - - ^„ 

will enable one to obtain direct information about the sub-leading contributions to the 
moments and, as a result, directly assess the quality of the equations for the moments 
Sn{r). The knowledge of /3„ will define universality classes, differing by the mechanisms of 
turbulence generation. According to this work, the functions F2n+i should demonstrate the 
scaling behaviour all the way up to the integral scale L. If this is indeed so, the measurements 
are not too difficult. 

If turbulence is driven by the large-scale body force, the subleading correction to the 
Kolmogorov relation for 6*3 (r) is an analytic O(r^) function. This fiow can be realized in 
numerical experiments. In the real-life situations this kind of forcing rarely exist. The 
appearence of the 0(r^'^) non-analytic correction in the relation for S^lr), derived from the 
equation (22), can be easily explained in both cases of decaying and sheared turbulence 
considered above. The measurements in jets and wakes are usually taken at the distance 
X from the origin (nozzles, bodies etc) where turbulence is produced. Thus, as was stated 
above, the proper model is that of decaying turbulence at the time T — x/U after turbulence 
generation at i = 0, where U is the mean velocity at the crossection x. Then, assuming a 
close-to- self-similar decay the correction is: 

0{rS2t{r)) ~ rS2{r)Ft{T) ^ r^-' 
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where the function F{t) describes the time-dependence of 5*2 (r,t) in decaying turbulence. In 
case of the shear-generated turbulence the correction is [18]: 

0{r—S,{r))^r'-' 

We do not know how general this result is. According to present work, it cannot be universal. 
The theory makes a direct connection to Landau's remark about the role of the large-scale 
fluctuations of turbulence production. We cannot answer the most important question about 
universality of the exponents: for this we need detailed and quantitative information on the 
production terms in the equations of motion. However, even if the exponents are universal 
or belong to some broad universality classes, the probability density function of velocity 
difference is not universal since the amplitudes of the moments 5'„(r) depend on some of 
the features of the non-universal large-scale dynamics and the details of the single-point 
probability density. 

The assumption about gaussian single point PDF, responsible for the values of the am- 
plitudes of the even-order moments, was used here as a good approximation sufficient for 
demonstration of the basic features of the theory. It is clear that the non-universality of 
the symmetric part of the PDF Pe{U, r) is determined by deviations of the single-point PDF 
from the gaussian. It is interesting that, according to the present work, even neglecting 
non-universality of Pe{U, r), the asymmetric part Po{U, r), responsible for the odd-order mo- 
ments, is not universal. This is quite reasonable since the very existence of the flux, reflected 
in Po{U, r), is the result of the large-scale dynamics. 

The theory, presented here, is a departure from all previous field-theoretical attempts to 
develop an infra-red divergence-free turbulence theory, able to explain anomalous scaling of 
the moments of velocity difference. It has always been assumed that due to galilean invari- 
ance, the vertex corrections are equal to zero in the infra-red limit k ^ 0. The supposed 
GI led to formulation of the Ward identities which were not too helpful. The low- order 
Kraichnan's LHDIA [28] , which in addition to the conservation laws correctly accounted for 
such basic symmetries of the problem as random galilean invariance, led to the Kolmogorov 
energy spectrum without any corrections. It is interesting, that the same approximation. 
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applied to Burgers turbulence, resulted in the k^^- energy spectrum corresponding to strong 
shocks. Kraichnan explained this non-trivial result in terms of the phase-decorrelation due 
to the interaction between components of the velocity field, non-existent in the Burgers 
dynamics, which effectively prevents the shock formation. According to [28] , it is this decor- 
relation which is responsible for the formation of the close- to- experimental data Kolmogorov 
5/3-energy spectrum. All attempts to preserve galilean invariance of the theory of three- 
dimensional turbulence led to disappearence of the integral scale L from the problem and re- 
sulting inability of the theory to predict deviations from the Kolmogorov scaling. Polyakov's 
[11], Boldyrev [29] and Parisi et. al.[16] theories of Burgers turbulence showed that the 
galilean invariance is not to be taken for granted: only the low-order moments Sn<i oc r" 
corresponding to \U\ << Urms can be described in this regime. The structure functions 
Sn oc r with n > 1 scale with Urms depending on the large- scale features of the flow. These 
works stated that galilean invariance is not sacred and departures from it are responsible for 
the anomalous scaling of the high-order moments observed in Burgers turbulence. The same 
conclusion was derived in an earlier work [30] predicting 

{u{x + r) - u{x))S{x)£{x + r) oc Urmsi^Yir / Lf 

This experimentally confirmed relation [31] explicitly involves the GI breaking Urms- 

The present work makes an additional step in this direction: it assumes that due to 
incompressibility, GI in three-dimensional turbulence is broken for all, even small velocity 
fluctuations. This means that the "normal scahng" does not hold for all moments n > —1 
which seems to agree with experimental data showing deviations from Kolmogorov scaling 
of all moments Sn- Free from the GI restrictions, the vertex corrections are introduced in 
this paper from the very beginning resulting in the equation of motion for the probability 
density of velocity difference. This equation is based on some general symmetry properties 
of the system and satisfies all known realizability constraints. 

The theory developed here is based on a few assumptions. First of all, it assumes the 
existence of a closed equation for longitudinal structure functions Snij-). This is a mere 
generalization of the Kolmogorov result for a particular case with n = 3. Physical grounds 
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for the equation for Sn are based on the fact the Navier- Stokes non-hnearities tend to produce 
two main effects: the shock-generation due to the advection terms which are balanced by 
the pressure contributions. An interplay between the two leads to creation of the vortical 
structures seen in the experimental data. The three-dimensional nature of the structures is 
lost when one considers projection of the entire dynamics onto a line. All we know is that 
the shock-production is effectively prevented by the pressure terms leading to invalidation 
of the bi-fractal description of the pure-Burgers dynamics. We also know that the equation 
of motions are invariant under transformation [/—>—[/ and x — > —x and that the dynamic 
equation for the A^-point generation function must satisfy the general fusion rule transforming 
it, upon point merging, into the equation for the N — 1-function. These are the physical 
reasons, responsible for all, but one, contributions to the equation (20). 

As was mentioned above, except for the A-teicm in (20), all others are more or less 
prescribed by the original equation of motion. However, without the 74-term, the equation 
(20) contradicts exact relation (42) and thus cannot be correct. I have not been able to find 
an alternative description, obeying general the fusion rules and symmetries of the problem, 
and producing solution satisfying (42). For example, one can add the 0(^2) -contribution to 
the right side of (19) which does not contradict the basic symmetries. It violates, however, 
one of the principle constraints of the theory 5"! =< U >— and thus, cannot be correct. 
The success of the Boldyrev theory including this term in description of some of the regimes 
of Burgers turbulence [28] is based on the fact that, as Polyakov's work, it describes only the 
moments Sn with n < 1 and the result Si = can be achieved using contributions coming 
from the non-scale invariant terms, which are beyond approximations of Refs.[ll].[29]. The 
same happens in the Parisi et. al theory [16] leading, in accord with the bi-fractal picture, 
to the PDF consisting of two contributions: responsible for the moments with n < 1 and 
n > 1, respectively. In the present paper we, treating all moments S'„(r) with n > on 
equal footing, do not have the luxury of satisfying the dynamical constraints using some 
contributions extraneous to the theory and, as a result, our choice of the allowed terms in 
the equations of motion is much more narrow. One may also attempt to add 
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not violating the symmetry of the equation (20). This gives 

_ n + /i 
n + B 

However, since = 0, the constant h — 0. 

The equation (22) for the PDF can be rewritten in the hmit of small r as: 

dP ^dP ^„ fU dP(y,r) , P 

where r oc r/U. The meaning of the A- term can be understood from the shape of (48) 
if we take into consideration that, unlike in the one- dimensional case, the interaction with 
the transverse components of the velocity field, produces an effective sources and sinks or 
friction for the longditudinal correlations. Then, (48) is the equation of motion taking 
these sinks and sources into account in the relaxation- time approximation. In other words 
T ^ r/U is simply the life-time (eddy turn-over time) taking the longditudinal structure to 
substantially change its shape and size due to interactions with the transverse components. 
This characteristic time, though plausible, is yet to be derived from the final "microscopic 
theory". It is important that the equation (48) is conservative, so that, since < U >— 0, 
So — 1 . Thus, the right side of (48) describes both sinks and sources. This is consistent with 
the results of the present paper: in the small-scale limit: S'„ oc r^" >> rt ^ SnB ^oi n < 3 
while 5'„ << SnB oc r for all n > 3. Here Sns is the n^'^-order moment of velocity difference 
measured in the large-scale-driven Burgers turbulence. This can be seen directly from the 
equation (48): the PDF tail with U > grows while the part those with U < decreases 
making the Navier-Stokes PDF P{U, r) much more symmetric than the one governing the 
Burgers dynamics. This means that due to the interaction between different components of 
the velocity field, the weak structures {\U\ < Urms), generated by the 3D the Navier-Stokes 
dynamics are much more intense than their counterparts in the Burgers turbulence. At the 
same time, due to Kraichnan's phase decorrelation, the strong structures in the Navier-Stokes 
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turbulence are much weaker than the strong shocks, responsible for the high-order moments 
of the Burgers dynamics. 

In a recent paper by Zikanov et.al. [32] a modified one-dimensional Burgers equation 



has been considered. As in Ref. [12], the white-in-time gaussian random force characterized 
by the spectrum \f{k)\^ oc was used. The large scale dissipation was introduced to avoid 
growth of the mode u{k = 0). It has been shown that addition of the non-local contribution 
is sufficient to prevent the shock formation and generation of the non-trivial exponents 7^ 1 
for n > 3. The possible relation of (49) to the equations, introduced in this paper, will be 
discussed elsewhere. 

The good agreement of the scaling exponents of the structure functions wiht experimental 
data, derived in this paper, though gratifying, is not the most important outcome. The 
cascade models, not related to the equations of motion, also give quantitatively correct ^2n- 
However, no model was able to address the problem of the asymmetry of the probability 
density function P{U,r) ^ P{—U,r) and, as a consequence, predict the scaling exponents 
and the amplitudes of the odd-order structure functions. Only dynamic theory, based on 
the Navier-Stokes equations, invariant under transformation u —u and x —x can lead 
to the asymmetric probability density and correct properties of the odd-order moments. 

The present theory, though based on some physical considerations developed for the 
Burgers equation, is not a small perturbation around Burgers phenomenology: the coeffi- 
cient B and B » h 2 obtained in [11]. Moreover, the relevant coefficient 
renormalizing advection contributions in equation (19), is strongly negative unlike parame- 
ter 6 > in the theory of Burgers turbulence. This means that the pressure and transverse 
terms, preventing formation of strong shocks, are extremely important here. On the other 
hand, the Burgers effects are not unimportant at all. To demonstrate it, we can neglect the 
"original" Burgers terms in the equation of motion and derive the relation for the moments: 




(49) 



34 



A „Kn_l+l 

Sn = ^nr- - ^^n{n - 1)^ 

where now 



n ^ B' 



_ A 

'^^ - ;g7^ 

and A/B' < 1. We see that without "small" Burgers terms the equation of motion gives 
normal scaling and that is why they are essential for derivation of the anomalous scaling 
exponents. If all this is correct, one may say that the anomalous scaling is the result of 
the dynamic interplay of the Burgers-like tendency to create singularities (shocks) with the 
"normalizing" action of the pressure terms and the incompressibility constraints. 

The equations developed here are based on the phenomenology relevant for the lon- 
gitudinal structure functions and that is why wc cannot say anything about shape and 
scaling of transverse structure functions. The main problem is that the odd-order mo- 
ments 5'2„+i =< (va(x-l-r) — Vo,(x))^"+-'^ >= 0, where vector r is parallel and is a 
component of the velocity field perpendicular to the x-axis. This means that the PDF 
P{Ava,r) = P{—Ava,r) and the equations governing probability density of transverse ve- 
locity differences must have different symmetry properties than (20)-(21). 

The most important feature of hydrodynamic turbulence, distinguishing it from equi- 
librium statistical mechanics, is constant energy flux in the wave-number space. It is this 
energy flux that makes the probability density P{U,r) asymmetric, leading to the non-zero 
values of the odd-order longitudinal structure functions. It is not clear how the information 
about the energy flux is reflected in the transverse structure functions coming out from the 
corresponding symmetric probability density. It is even unclear if S^{r) are a dynamically 
relevant object. In the theory presented here, the transverse components of the velocity 
fleld simply serve as a "bath" introducing some renormalization and dephasing into the the 
energy flux- carrying longitudinal dynamics. The angular dependence of the structure func- 
tions in three-dimensional turbulence can be recovered using the multidimensional equation 
for the two-point generating function with the pressure terms accounted for in the mean fleld 
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approximation, similar to the one introduced in this paper. It is not clear if this can lead to 
the improved description of experimental data. 

Acknowledgments I am grateful to R.H. Kraichnan whose remarks, insights and con- 
structive suggestions were most essential for this work. My thanks are due to S. -Y.Chen, 
K.R.Sreenivasan, B. Dhruva, R.Miles and U.Frisch for providing me with their most recent, 
sometimes unpublished, experimental data. Very interesting and stimulating discussions 
with S. Boldyrev, A. Chekhlov, M. Chertkov, U. Frisch, M. Nelkin, A.Polyakov and B. 
Shraiman are gratefully acknowledged. This work was supported in part by the ONR/URI 
grants, 
references 

1. L.D. Landau and E.M. Lifshitz, Fluid Mechanics, Pergamon Press, Oxford, 1987; A.S.Monin 
and A.M.Yaglom, "Statistical Fluid Mechanics" vol. 1, MIT Press, Cambridge, MA (1971) 

2. A.N. Kolmogorov, J. Fluid Mech., 5, 497 (1962). An interesting chain of events starting 
with Landau's remark and ending with the Kolmogorov paper is described in: U. Frisch, 
"Turbulence", Cambridge University Press, 1995. 

3. R.H. Kraichnan, Phys.Rev.Lett.,72, 1016 (1994) 

4. F. Gawedzki and A. Kupianen, Phy.Rev.Lett., 75, 3834 (1995) 

5. M. Chertkov, G. Falkovich, I. Kolokolov and V.Lebedev, Phys.Rev. E,52, 4924 (1995) 

6. M. Chertkov, G. Falkovich, I. Kolokolov and V.Lebedev, Phys.Rev. E,51, 5609 (1995). 

7. B. Shraiman and E. Siggia, CR Acad.Sci., 321, Serie II (1995) 

8. V.Yakhot, Phys. Rev. E, (1996) 

9. M.Vergassola and A.Mazzino, |chao-dyn/9702014 

10. O.Gat, V. L'vov, E. Podivilov and I. Procaccia |chao-dyn/9610016 

11. A.M. Polyakov, Phys.Rev. E 52, 6183 (1995) 

12. A. Chekhlov and V. Yakhot, Phys.Rev.E 51, R2739 (1995) 

13. A. Chekhlov and V. Yakhot, Phys.Rev. E 52, 5681 (1995) 

14. V. Yakhot and A.Chekhlov, Phys.Rev.Lett. 77, 3118 (1996) 

15. M. Chertkov, I. Kolokolov and M. Vergassola, preprint 

16. J.P. Bouchaud, M. Mezard and G. Parisi, Phys.Rev. E 52 3656 (1995); V.Gurarie 



36 



and A.Migdal, Phys. Rev. E54,4908 (1996); the relation for Si{r,9) has been written by 
M.Nelkin; A.Chekhlov and V.Yakhot, unpubhshed (1996) 

17. D. Forster, D. Nelson and M. Stephen, Phys. Rev. A16 732 (1977) 

18. S. Corrsin, NACA R&M 58B1 (1958); J.L. Lumley, Phys. Fluids 10, 855, (1967); I am 
grateful to B. Shraiman who brought my attention to these works. 

19. S.-Y. Chen, private communication, 1997 

20. M. Chertkov, Phys.Rev. E, 55, 2722 (1997) 

21. P. Kailasnath, A.A. Migdal, K.R. Sreenivasan, V. Yakhot and L. Zubair, Yale preprint 
(1991) 

22. S.-Y. Chen, private communication, 1997 

23. B. Dhruva and K.R. Sreenivasan, private communication (1997) 

24. A. NouUez, G. Wallace, W. Lempert, R.B. Miles and U.Frisch, J.Fluid.Mech. 339, 287 
(1997) 

25. Ya. G. Sinai and V.Yakhot, Phys.Rev.Lett., 63,1962 (1989) 

26. S. Pope, Combust. Flame 27, 299 (1976) 

27. I am gratefull to R.H.Kraichnan who pointed out this relation to me 

28. R. H. Kraichnan, Phys.Fluids. 11, 266 (1968) 

29. S. Boldyrev, Phys.Rev.E 55, 6907 (1997) 

30. V. Yakhot, Phys.rev. E 50, 20 (1994) 

31. A. Praskovskii and S. Oncley, Phys. Rev. E 51, R5197 (1995) 

32. O. Zikanov, A. Thess and R. Grauer, Phys.Fluids 9, 1362 (1997) 



37 



Figure legends 

Fig. 1. Comparison of the calculated scaling exponents (formula (27)) with the results of 
numerical simulations of Chen [19]. 



Fig. 3. Probability density P{U,r) as a function of U/urms- 

Fig. 4. Mcsured PDF's of transverse structure functions [24] for a few values of the displace- 
ment r = 3600/xm; 900/xm; 28/im;. 

Fig. 5. Approximate parametrization of the PDF r°'P{U,r) ys x — ^ using formula (37). 
Fig. 6. Asymmetric PDF P{U,r) given by (36)-(37). 

Fig. 7. Calculated conditional mean ^^^''"^ using approximate expression for the PDF (36)- 



Fig. 2. Probabihty density of velocity differences F{x) — rsP^^-) vs x 
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Figure 1: Comparison of the calculated scaling exponents (formula (27)) with the results of 
numerical simulations of Chen [19] 
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Figure 2: Probability density of velocity differences F{x) — r3P(-^) ys x = ^- 

r3 r3 
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Figure 3: Probability density P(U,r) as a function of U 
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Figure 5: Approximate parametrization of the PDF r"P{U, r) as a vs a; = ^ using formula 
(37) 
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Figure 6: Asymmetric PDF P(U,r) given by (31)-(32) 
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Figure 7: Calculated conditional mean ^^'^■^ using approximate expression for the PDF 
(36)-(37) 
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Figure 8: Calculated normalized moments Ss{r) and S'5(r) 
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